ON THE GLOBAL REGULARITY OF 
SUB-CRITICAL EULER-POISSON EQUATIONS WITH PRESSURE 



EITAN TADMOR AND DONGMING WEI 

Abstract. We prove that the one-dimensional Euler-Poisson system driven by the Poisson forcing 
together with the usual 7-law pressure, 7 > 1, admits global solutions for a large class of initial data. 
Thus, the Poisson forcing regularizes the generic finite-time breakdown in the 2x2 p-system. Global 
regularity is shown to depend on whether or not the initial configuration of the Riemann invariants 
and density crosses an intrinsic critical threshold. 



1. Introduction 

It is well known that the systems of Euler equations for compressible flows can and will breakdown 
at a finite time even if the initial data are smooth. A prototype example for such systems is provided 
by the 2x2 system of isentropic gas dynamics 



(1.1) 



pt + (pu) x = 
(pu) t + (pu 2 ) x = -p x 



where the pressure p = p{p) is given by the usual 7-law, p{p) = Ap 1 . By using the method introduced 
in [La64j to deal with pairs of conservation laws, it can be shown that (jl.ip will lose the C 1 -smoothness 
due to the appearance of shock discontinuities unless its two Riemann invariants are nondecr easing. 
Thus, the finite time breakdown of (11. ip is generic in the sense that it holds for all but a "small set" 
of initial data. 

On the other hand, if we replace the pressure by Poisson forcing, then we arrive at the system of 
Euler-Poisson equations 

- ^ I Pt + (pu)x = 0, 

[ {pu) t + {pu 2 ) x = -kptp x k > , 

subject to initial data (no , po > 0) . Here ip = (p(p) is the potential, which is dictated by the (one- 
dimensional) Poisson equation, ip xx = — p . In this case, there is a "large set" of initial configurations 
which yield global smooth solutions. More precisely, [ELT01| have shown that (|1.2fl admits a global 
smooth solution if and only if 



(1.3) u 0x (x) > -y/2kp (x) . 

Thus, following the terminology of [LT02), the curve uq x + ^/2/cpo = is a "critical threshold" in 
configuration space which separates between initial configurations leading to finite time breakdown 
and a "large set" of sub-critical initial configurations which yield global smooth solutions. In particular, 
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(|l,3p allows negative velocity gradients (- depending on the local amplitude of the density), which 
otherwise are excluded in the case of inviscid Burgers equations, corresponding to k = 0. 

In this paper we turn our attention to the full Euler-Poisson equations driven by both - pressure 
and Poisson forcing, 

IPt + (pu) x = , 
{pu) t + {pu 2 ) x = -p(p) x - kptfx, k>0, 
-<Pxx = P ■ 

These equations govern different phenomena, ranging from the largest scale of e.g., the evolution 
gravitational collapse in stars, to applications in the smallest scale of e.g., semi-conductors. There is 
a considerable amount of literature available on the local and global behavior of Euler-Poisson and 
related problems. Consult [Ma86] for local existence in the small .£P-neighborhood of a steady state 
of self-gravitating stars, |CW96] for global existence of weak solutions with geometrical symmetry, 
[Gu98] for global existence for 3-D irrotational flow, [MN95J for isentropic case, and |JR00] [PRV95J 
for isothermal case. Consult [Pe90] |MP90j [Sigg] [En96] |WC98| . |BW98] and in particular, [En96] 
(more about that below), for non-existence results and singularity formation. The question of global 
smoothness vs. finite breakdown was studied in a recent series of works of Engelberg, Liu and Tadmor, 
in terms of a critical threshold phenomena for 1-D "pressure-less" Euler-Poisson equations, [E LTOlj 
and 2-D restricted Euler-Poisson equations, 1.T021 |LT03| . 

The natural question that arises in the present context of full Euler-Poisson equations (jl.4p is 
whether the pressure enforces a generic finite time breakdown or, whether the presence of Poisson 
forcing preserves global regularity for a "large" set of initial configurations. We answer this question 
of "competition" between pressure and Poisson forcing, proving that the Euler-Poisson equations (jl.4p 
with 7 > 1 admit global smooth solutions for a "large set" of sub-critical initial data such that 

(1.5) u 0x (x) > -ETo7^M+ v^7 IPOx( ^i , 7>1- 

po(x) 2 

Here, Kq is a constant depending on k, 7 and the initial data. In the particular (and important) case 
of isothermal equations, 7 = 1, we have Kq = y/2k and (|l,5p amounts to a sharp critical threshold, 

(1.6) n ,(x)>- V / 2fc^ + v^^M, 7=1- 

Po{x) 

The inequalities (|1.5p . (|1.6p quantify the competition between the destabilizing pressure effects, as the 
range of sub-critical initial configurations shrinks with the growth of the amplitude of the pressure, 
A, while the stabilizing effect of the Poisson forcing increases the sub-critical range with a growing k. 
In particular, (|1.6p with A = recovers the pressure- free critical threshold (jl.3p . 

Formation of singularities and global regularity of (]1.4p were addressed earlier by Engelberg in 

[En96]. His results show finite-time break-down if uq x (x) — \/A~f\po x (x)\po(x) is "...sufficiently 
negative at some point". Our contribution here is to quantify the critical threshold behind this 
asymptotic statement. To fully appreciate this quantified threshold, we turn to the converse statement 

in [En961 Theorem 2]: it asserts the global regularity of (jl.4p for a class of initial data such that 

7-3 

uo x (x) — T/Ay\po x (x)\po(x) 2 > 0. It is a "non-generic" class (in the sense of requiring both Riemann 
invariants at t = to be globally increasing). In fact, by (|1.5p one has a negative threshold, — Ko^/po, 
implying the existence of a "large" class of sub-critical initial data with global regularity. 

The paper is organized as follows. In section [21 we reformulate the system f|l .4[) with its Riemann 
invariants as a preparation for the analysis carried out in sections [3] and [H In section [31 we prove 



ON THE GLOBAL REGULARITY OF SUB-CRITICAL EULER-POISSON EQUATIONS 3 

our main results, providing sufficient conditions for "large sets" of sub-critical initial configurations 
which yield global smooth solution. In section HI we give examples of finite time breakdown for super- 
critical initial data. Combining our results in sections [3] and HI they confirm the existence of a critical 
threshold phenomena for the full Euler-Poisson equations (jl.4p . 

2. Riemann Invariants 

2.1. The Euler-Poisson equations with 7-law pressure: 7 > 1. We begin by rewriting the 
Euler-Poisson equations (jl.4p as a first order quasilinear system 

(:).*' 

u p 

where the Jacobian J := | | has two different eigenvalues 

A-fp^- 2 u 





A := u — a/ A^p~^~ < p := u + y A^p~^~ . 
and let R and S denote the Riemann invariants of the corresponding Euler system (jl.ip 

(2.2) R:=u- 2 -^p^ and S := u + ^ p 2 ^ . 

7 — 1 7 — 1 

They satisfy the coupled system of equations, 

(2.3a) Rt + XR X = -kip x , 

(2.3b) St + pS x = —kipx , 

coupled through the Poisson equation — 4> xx = p . If we set r := R x , s := S x then upon differentiation 
of (g3D we get 

(2.4a) r t + \r x + \ s rs + \ R r 2 = kp , 

(2.4b) s t + ps x + p s s 2 + Prts = kp . 

Next, we observe that A = — ^^^-(S — R) and p = + : ^(<S — R). Hence, expressed in terms 
7 - 1 

of 8 := — - — , we have for 7 > 1, 

1 + 61 , , 1-0 .7-1 

*R = PS = — Y~ and X S = PR= — Y~ ' 6 := - °' 

and the pair of equations (j2.4j) is recast into the form 

/ _ y. , 1+9 2 1-6 

(2.5a) r H — r H — rs = fcp, 

, ^ , l + 6» 2 1-0 

(2.5b) s + -^-s l + ~^-rs = kp . 

Here and below {} v := dt + Xd x and {}' := dt + pd x denote differentiation along the A and p particle 
paths, 

T A := {(x,t) I x(t) = X(p(x,t),u(x,t))}, T M := {(x,t) \ x(t) = p(p(x,t),u(x,t))} . 
To continue, we rewrite the equation for p as 

(2-6) (pt + Xp x ) + ^-^p x + p^-^ = 0. 
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Since s — r = S x — R x = 2\J A'yp""^ p x ; it enables us to express ^— — p. 



7-1 s — r 

Ajp 2 Px = p—^- , so 



2 v 2 

that the equation (12. 6h can be written along the A particle path as p K + ps = 0. Similarly, it can be 

written along the p particle path as p' + pr = 0. Assembling the above equations together, we arrive 

at the following system governing r, s and p, 



(2.7a) 

and 

(2.7b) 



1 + 



-r z + 



1 



-rs 



kp, 



p' + ps = , 



' / 1+0 2 1 

s + — — s + - 



-rs 



2 

p' + pr 



kp , 
0. 



Finally, we use the integration factors 1/y/p and r/2p^/p in the first and second equations of each pair 
in (12.71). to conclude 



ky/p, 





) + 2 


J r 2 


6 rs 


s/Pj 




2^ 


s ^ 


\' 1 + ^ 
) + 2 


J s 2 


rs 






2^ 



(2.8a) 
(2.8b) 



2.2. The isothermal case 7=1. In this case, the two eigenvalues are A = u — \f~A < p = u + \[A 
with the corresponding Riemann invariants R = u — \f~Ah\ p and S = u + \AAln p. Their derivatives, 
r and s, satisfy the pair of equations, corresponding to (|2.8aj) . (|2.8b[) with 9 = (7 — l)/2 = 0, 



(2.9a) 
(2.9b) 



( 



r A 1 r 2 

(^) + 5^"*^ 



3. Global smooth solutions for sub critical initial data 

For the pressure-less Euler-Poisson equations (|1.2p . the evolution of u x and /? could be traced 
backwards along the same particle path to their initial data at t = 0. The scenario becomes more 
complicated with the additional pressure term, due to the coupling of r and s along different particle 
paths which are traced back to different neighborhoods of the initial line t = 0. This is the main 
obstacle in finding the sharp critical threshold of the full Euler-Poisson system (|1.4p . To this end, we 
will seek invariant regions for the coupled system, governing the Riemann invariants. We begin this 
section with the following lemma. 

Lemma 3.1. Given that the total charge Eq := po(x)dx < 00 , then p(x,t) and u(x,t) remain 
uniformly bounded for all t > 0. 

Proof. Under the given condition, we can set (e.g., [ELT014 p. 116]) 
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which satisfies — Eq < (p x (x, t) < Eq , for all t > and x £ R . 

Recall the transport equations (|2.3ap . (l2.3bjl which govern the Riemann invariants along different char- 
acteristics R}-\-ktp x = S' + kip x = 0. Since (p x is bounded, these transport equations tell us that R and 
S remain uniformly bounded with at most a linear growth in time. Indeed, for all M 3> 1 we have 

(3.1) sup {\R(x,t)\,\S(x,t)\} <C + kE Q t, C := sup {\R (x)\,\S (x)\} . 

\x\<M ^^M+Uoot 

Take the sum and difference of S and R to find that u(x, t) and p(x, t) in (|2.2p remain bounded, 

{(C + kE t)^, 7>1, 
exp(kE t), 7=1. 

□ 

We note in passing that the time growth asserted in (|3.2p is probably not sharp; the estimate can be 
improved after taking into account the uniform bounds of R x j ^fp and S x / ' yfp discussed in theorems 
O and O below. 

Remark 3.1. According to Lemma \3.1l the only way that the full Euler-Poisson system fil.4\ ) breaks 
down at a finite time is through the formation of shock discontinuities where \u x \ and/or \p x \ blow up 
| oo , but neither will concentrate at any critical point. This is in contrast to the breakdown of the 
"pressure-less" Euler-Poisson equations fl.ty) . where —u x (x,t) = p(x,t) f oo simultaneously at the 
critical time. 

3.1. Critical Threshold for isothermal case: 7=1. We begin with the isothermal case, 7 = 1, 
which plays an important role in various applications. Compared with the general case (|2.8p . the 
isothermal case becomes simpler due to the fact that 8 = decouples the dependence on r and 
s through the mixed term 9rs, which disappears from left hand side of (|2.9p . Here we prove the 
following sharp characterization of the critical threshold phenomena. 

Theorem 3.1. Consider the isothermal Euler-Poisson system |j.^[ ) with pressure forcing p(p) = Ap , 
and subject to initial data (uq,pq > 0) with finite total charge, Eq = po(x)dx < 00. The system 
admits a global smooth, C 1 -solution if and only if 

(3.3) u 0x {x) > -v / 2fcp (x) + y/A ^fr > VxER. 

Remark 3.2. Expressed in terms of the Riemann invariants specified in j pTj[ u x ± \fAp x j p, theorem 
\3.1\ states that the isothermal Euler-Poisson equations admit global smooth solutions for sub-critical 
initial conditions, 



(3.4) so > —\/2kpo and ro > —yZkpQ. 

Proof. We define X := — and Y := . Equations (I2.9aj) . ()2.9bh then read 

sJ~P V~P 

(3.5a) X x = ^-{2k - X 2 ), 

(3.5b) Y' = ^-{2k-Y 2 ). 
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It follows that 



xH 



> 0, X £ 

= 0, IXI = 



< 0. 



2k,V2k] 
2k, 
X\ > V2k, 



and similarly, 



> o, y e 

=0, \Y\: 



< 0. 



2k,V2k 
2k, 
Y\ > V2k. 



Thus, starting with (|3.4|) . Xq, Yq > — \[2~k, we find that X and Y remain bounded within the invariant 
region [—y/2k, V2k], or otherwise, they are decreasing outside this interval. We conclude that 



X(-,t),Y(;t) < max{V2fc,X o (-)To(0} 



Lemma 13.11 tells us that p is bounded. The boundedness of X , Y and p imply that r = X^fp and 
s = Y \J~p~ remain bounded for all t < oo, and hence the Euler-Poisson system (jl.4p admits a global 
smooth C 1 -solution. 

Conversely, suppose that there exists Xq = X(xq) < —\[2~k. We will show that this value will evolve 
along r\(xo,0) such that X(-,t) will tend to — oo at a finite time. To this end, assume that Y is well 

behaved, i.e., Yq(-) > —y/2k so that Y(-,t) < Y% := max |Yo(-), \/2fc| for all t's (otherwise, the finite 
time blow up of Y can be argued along the same lines). It follows that s = Y \J~p~ < Y\J~p and inserting 
this into p s = —ps, we find p v > —Y\p^l 2 . This yields the lower-bound 



P> 



Yit + 2/^po- 



and together with (|3.5aj) . we conclude that X{-,t) satisfies the following Ricatti equation along the 
IVpath, 



(3.6) 



Integration of (|3.6p yields 
(3.7) 

Thus, starting with Xq < 
X(t 1 t c ) tends to — oo. 



Fit + 



X(;t) < 



X' 



X 1 := (Xl - 2k) /Xl > 0. 



Yt 



Xi In (1 + y/p^t/2) + YiX 

2k < it follows that there exists a finite critical time t c > such that 

□ 



The critical threshold condition (|3.3p reflects the competition between the Poisson forcing and the 
pressure. It yields global smooth solutions for a "large" set of initial configurations allowing negative 
velocity gradients. In the particular case that there is no pressure, ^4 = 0, (|3.3p is reduced to the critical 
threshold condition of the "pressure-less" Euler-Poisson equations uq x > — y / 2kpo(x) of (ELTOlj . 
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3.2. Critical threshold for 7 > 1. The equations for the Riemann invariants (|2,8ap . (|2.8bp are 
coupled through the mixed term, Ors/2. We note in passing that it is possible to get rid of this mixed 
term when integrating (|2.7aj) . (|2.7b|) with the integration factors p^ -3 )/ 4 , and rp^~ 7 ^ i (3 — j)/4 in 
the first and second equations in each pair, yielding 



9-1 \ 1 + 9 o 1+2 
rp 2 J H — r p 2 = kp 2 , 

e-i \ 1 + 9 o tzl 1+1 
sp 2 I -| - — s p 2 = kp 2 . 

Nevertheless, it will prove useful to use the same integration factors, 1/^/p and r/2p^J~p which led 
to (|2,8p . The main task is to identify the invariant region associated with (|2.8p , corresponding to the 
isothermal invariant region [—y/2k,y/2k] discussed in theorem 13.11 

Theorem 3.2. Consider the Euler-Poisson system with 7 law pressure p(p) = Ap 1 , 7 > 1, 

subject to initial data (uo,po > 0) with finite total charge, Eq = J°° po(x)dx < 00. Then, there exists 
a constant Kq > depending on fe,7 and the initial conditions (specified in 113. 9b\) below), such that 
the Euler-Poisson equations admit a global smooth, C 1 -solution if, 

I pox 0)1 



(3.8) u Qx (x) > -K y/p (x) + y/Ay- . 

p (x) 2 

Before we turn to the proof of this theorem, several remarks are in order. 

Remark 3.3. Expressed in terms of the Riemann invariants, r = u x — a/ ' A^p§ x j p^ 7 ^ 2 and s = 
u x + \/A^ypo x /p^ 7 ^ 2 ; the critical threshold i3. 8\) reads 

( o n \ r (a;) s (x) 

(3.9a) , — > — Kp. 

s/po(x) \/po(x) 

The constant Kq is given by 

(3.9b) K = u , M = max <^ V2k, r == , r == > 

2(1 + 0) x I s/po(x) Vpo( x ) J 

We mention two simplifications which are summarized in the following two corollaries. We first 
observe that if the initial configurations satisfy the upper-bound ro(x), sq(x) < ^j2kpo(x) then (|3.9bp 



y2k 

yields Mq = V2k, hence Kq = -, and theorem 13.21 implies the following. 

1 + 9 



Corollary 3.1. Consider the Euler-Poisson system |j.^[ ) with 7 law pressure p(p) = Ap 1 , 7 > 1, 
subject to initial data (uq,pq > 0) with finite total charge, Eq = po(x)dx < 00. Then, the Euler- 
Poisson equations admit a global smooth, C 1 -solution if for all x £ R, 

(3.10) |«oxOc)| < v / 2kp~ Jxj - ^ lPOx{x ^ . 

p (x) 2 

9 M - (9 M + J8k(l + 9)) /y/2 

The next result follows from the trivial inequality — Ko < ; . 

M J u ~ 2(1 + 6) 

Corollary 3.2. Consider the Euler-Poisson system |j.^[ ) with a ^y-law pressure p(p) = Ap' 7 , 7 > 1, 
subject to initial data (uo,po > 0) with finite total charge, Eq = J po(x)dx < 00. Then, the Euler- 
Poisson equations admit a global smooth, C 1 -solution, if for all x G R, 
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V^/ 2(7 + 1) a [ po^-^J p x 



max 



3-7 
2 



Remark 3.4. We observe that as in the isothermal case, the critical threshold in its various versions 
\3. 8\) , nnfy . 113. 10]) and 113.11]) , allow a "large set" of initial configurations with negative velocity gradi- 
ent, due to the competition between the stabilizing Poisson forcing kp<p{p) x and the destabilizing pres- 
sure A{p~ l ) x . In the extreme case that Poisson forcing is missing k = , the breakdown of the system 
is generic unless uq x is positive enough (so that ro, so > 0). In the other extreme of a "pressure-less" 
Euler-P oisson, A = 0,7 = 1, the critical thresholds i3.8\) . $3.10]) are reduced to uq x {x) > — ^/2kpo(x) , 
which coincides with the "pressure-less" critical threshold hi. 3]) found in [ELTOlj . 

T S 

Proof. Expressed in terms of X := and Y := — , equations (|2.8p read 

Vp Vp 

(3.12a) X" = Jp(k - l -^~X 2 + °-XY 

(3.12b) Y' = J- P (k - i±V + °-XY 

We seek an invariant region of the form [— Kq, Mn], with Kq, Mq > yet to be determined. To this 
end we construct a "buffer zone" in which positive values of X, Y must decrease and hence remain 
upper-bounded. We begin by noticing that if X, Y < M thenS X + Y < M 2 , and recalling that 6 > 0, 
(ptZLZl) then yields 

(3.13a) X s < y/pfk - ^Y~ x2 + \ M ^j > X > °> 

(3.13b) Y '<^p[k- l -±iY 2 + e -M^y Y>0. 

This in turn implies that 

X or Y are decreasing whenever X G Xu or Y 6 Im- 

Here X M is the interval T M := (C+(M),M) where C+(M) := v / (2k + 6M 2 )/(1 + 9) is dictated by 
the largest root of the quadratics on the right of (|3.13j) . Observe that for Xm to be nonempty requires 



M > V2k. We therefore set, M e := max, j\/2fc+e, X (x), Y (x) >. We claim that X,Y < M e : indeed, 

either X, Y < C + (M e ) < M e or, if X, Y > C + (M e ), then they must decrease being "traped" inside 
Im e and hence X, Y < M e . Letting e { we end up with the upper-bound 

(3.14) X(-,t),Y(;t)<M Q , M :=m^W2k,X {x),Y {x)\. 



We let Z+ = max{X, 0} and Z- — min{Z, 0} denote the positive and negative part of Z. 
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In a similar manner, we study the lower bound of the invariant region. By (|3.14p and f)3. 12|) yield 
(3.15a) X" > Jp(k - ^-X 2 + 9 -M Q x\ , X < 0, 

(3.15b) Y' > Jp(k-l±iY 2 + 9 -M Q Y^ Y<0, 

which in turn, imply that 

(3.16a) X and Y are increasing if > X, Y > — Kq, 

where Kq is the smallest root of the quadratics on the right of f)3. 15j) , 



(3.16b) K := -OMo + ^Mi + mi+e) 

y ' u 2(1 + 9) 

The critical threshold condition (|3.8|) tells us that at t = 0, Xq,Yq > —Kq and (|3.16ap implies that 
X(-,t) and Y(-,t) remain above the same lower-bound, (|3.8p . As before, the bounds of X, Y and p 
imply that r = X^fp and s = Y^fp remain bounded, and hence the Euler-Poisson system (jl.4p a 
global smooth, C 1 -solution. □ 

4. Finite time breakdown for super-critical initial data 



Consider the Euler-Poisson system (II. 4p with a 7-law pressure, 7 > 1, and subject to initial data 
such that tq(x), sq(x) < \plk. Then, according to corollary 13 .1\ the following critical threshold is 
sufficient for the existence of global smooth solutions, 

Pq(x) 2 

In this section we show that this critical threshold is also necessary for global regularity. 



Theorem 4.1. Consider the Euler-Poisson system with a ^-law pressure p(p) = Ap 1 , 7 > 1, 

subject to initial data (uq,pq > 0) . The system loses the C 1 -smoothness if there exists an x £ M. such 
that 

\Pox{x)\ 



(4.1) uq x {x) < -^2kp (x) + v^47 



P0{X) 2 



Remark 4.1. Expressed in terms of the Riemann invariants, r = u x — y ' A^pQ X j p^ 7 ^ 2 and s 
u x + y/A^pox/pQ 3 7 ^ 2 ; the condition J^. 1\ ) reads 



(4.2) 3x£R s.t. r (x) < -y / 2k / 9 (x) , or s (x) < -y/2kpo(x) . 

The lack of smoothness in this case was shown in theorem I3.il for 7 = 1 and is extended for 7 > 1 
below. 

Proof. Recall equations (I3.12p for X : = and Y : = 

(4.3a) X^ = ^p(k- l -±lx 2 + 9 -XYy 

(4.3b) Y' = ^p(k- l -±^Y 2 + °-XYy 
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In the proof of theorem 13.21 we have shown that X and Y have an upper bound 
(4.4) X{-,t),Y(-,t) <M , M :=maxiV2k,X (x),Y (x)\. 



Suppose that there exists Xq = X(xq) < — \j2k. We will show that this value will evolve along 
I\(xo,0) such that X(-,t) will tend to — oo at a finite time. To this end, assume that Y is well 
behaved, i.e., 5^(0 > —v2~k so that Y(-,t) < M for all tfs (otherwise, the finite time blow up of Y 
can be argued along the same lines). It follows that along T\(xo,0) 

(4.5) r = Jp(k - i±^X 2 + 6 -XY^ < Jp(k - l -X 2 
Following exactly what we have done in the proof of theorem 13.11 we obtain the inequality 

(4.6) X(-,t)< 

where X% := (Xq — 2k) /Xq > . Thus, starting with Xq < —\f2~k < it follows that there exists a 
finite critical time t c > such that X(t f t c ) tends to — oo. □ 

We conclude with an example for a finite time breakdown. 
Example: Suppose at t = , uq(x) = and 



Xt In (1 + ^M t/2) + M X ' 



Po(x) 



1 , x < . 

x 

Ye 



1 - — , 0<x<e. 



l 

2 • 



x > e . 



Thus 

s fx) 



-,/ZFi{l-—)/2e, 0<x<e, 
, elsewhere . 



If we choose e small enough, then so (a;) < — y / 2kpo(x) for < x < e . According to theorem 14. 1\ the 
system (|1-4|) will break down in a finite time. This example shows that even if the fluid is near rest 
at t = 0, the pressure itself could still lead to collision. 
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